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f^ ■ We have measured the ac-Stark shift of the 4s ^5*0 — 4s4p Pi hue in Ca for perturbing laser 

f^ ' wavelengths between 780 nm and 1064 nm with a time domain Ramsey-Borde atom interferometer. 

^\l , We found a zero crossing of the shift for the ms = — mp = transition and a polarized perturbation 

. . . at 800.8(22) nm. The data was analyzed by a model deriving the energy shift from known transition 

" ' wavelengths and strengths. To fit our data, we adjusted the Einstein A coefficients of the 4s3d ^D — 
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4s4p P and 4s5s S — 4s4p P fine structure multiplets. With these we can predict vanishing 
ac-Stark shifts for the ^5*0 m — — ^Pi m = 1 transition and a~ light at 983(12) nm and at 
735.5(20) nm for the transition to the ^Po level. 

PACS numbers: 32.60.-|-i, 03.75.Dg, 32.70.Cs, 32.80.Pj 



I. INTRODUCTION 

o 

j^ ' The realization of the SI second by the ground state hyperfine transition in Cs yields the by far best realization 

^ of all SI units. Nevertheless, some optical frequency standards based on single ion traps [1, 2] or cold neutral atom 

Q . clouds [3, 4] are about to exceed the best microwave clocks in stability and accuracy. Due to the improved capabilities 

'^ ■ of frequency transfer with fs-combs [5, 6] these standards can be evaluated more easily. Thus fast progress can be 

J>^, expected. Recently, Katori and coworkers proposed and started to explore a so called "optical lattice clock" [7-9] . It 

^S^ combines the advantages of ion and neutral atom frequency standards by using large particle numbers and a confining 

Mh, optical potential, which allows long interrogation times and eliminates the Doppler effect and related uncertainty 

contributions. The wavelength and polarization of the optical trapping field are tailored such that both levels of 

the clock transition have the same Stark shift. Hence, the transition energy can be observed without perturbations 

^ . by the trapping field. The exact parameters have to be determined experimentally with high accuracy to meet the 

OO ' requirements of a lattice clock, since up to now calculations do not reach the necessary precision. 

In this article we present measurements of the Stark shift of the '*°Ca 4s^ ^Sq —isAp ^Pi transition to determine wave- 
lengths leading to vanishing frequency shift of the intercombination line. Our data is analyzed by fitting appropriate 
oscillator strengths to match measured and calculated Stark shifts. The transition probabilities determined in this 
^b ' way are compared to known values, which is especially of interest in cases where only theoretical results are available. 

o' 



II. THEORETICAL DESCRIPTION 



C/3 

^~j' The energy shift Ui{uj,p,mi) of an atomic state i with energy Ei and Zeeman level rui, which is induced by a 

tG . perturbing laser field with frequency v = uj/2Tr, polarization p, and irradiance / can be expressed in second-order 

^H' perturbation theory as Ui{u),p, rrii) — --ai{uj,p, mi)I/2eoc with the induced polarizability a^. We calculated the latter 
by summing the contributions from all dipole transitions from the desired state i to levels k with the respective Einstein 

k>( coefficients A^i (spontaneous emission rate for Ek > Ei), Zeeman levels m' , and transition frequencies z^i/j = LOik/^Tr 

^ [10] 

^ fi 3 ^ Ak^{2Jk + ll f J, 1 J, y 

The expression in large brackets denotes a 3J-symbol [11]. It describes the selection rules and relative strengths 
of the transitions depending on the involved angular momenta J, their projections m and the polarization p. The 
transition frequencies and Aki coefficients were taken from the comprehensive data collection of Kurucz [12]. Data 
from other sources, both, theoretical [15-18] and experimental [14, 19, 23] will be discussed in Sec. IV. Additionally, 
the influence of the continuum and highly excited states was approximated by using hydrogen wavefunctions [20] 
similar to the method of von Oppen [21]. 
A simplified level scheme is shown in Fig. 1. In Fig. 2 the function Ui{uj,p,mi) is depicted for different polarizations. 



Zeeman levels, selected levels and wavelengths. One observes that at the so called "magic wavelength" [9] (for 
appropriate polarization of the perturbing field) the shifts of 4s^ ^5*0 and 4s4p '^Pi states are the same and thus no 
frequency shift of the intercombination line is observed. This occurs for a polarized light and the mp = Zeeman 
sublevel of the '^Pi level around 800 nm, and at about 1000 nm for the mp = +1 (—1) level and a~ {(t~^) light. In 
the fermionic isotope, the transition '^Pq ~^ Sq is allowed due to hypcrfinc level mixing. For this transition a "magic 
wavelength" exists around 735 nm. 

The ac-Stark shift of the ■^Pi and ^Po levels is dominated in the depicted wavelength interval by the 1.9 /im, 612 nm, 
and 443 nm transitions shown in Fig. 1. Due to a vanishing value of the 3J-symbol in Eq. 1 some combinations 
of Zeeman levels and polarizations do not couple to the '^Si level by the 612 nm transition. E.g. in the case of 
(7~-polarization, the m = — 1 sublevel only couples to a to = —2 sublevel, which is not present in the '^5i level and 
thus the energy-shift Ui{u!,p,'mi) for this combination has no pole at that wavelength. 

For the prediction of the "magic wavelengths" , the 1.9 /im transition is of special interest, as its transition wavelength 
is the only one larger than the "magic wavelengths". Thus, only this line introduces negative contributions to 
ai{uj,p,mi) in Eq. 1. Therefore, its transition strength has, together with the 612 nm transition, a crucial influence 
on the slope of the curves in Fig. 2. The "magic wavelengths", depicted by circles in Fig. 2, are sensitively influenced 
by the counterplay of both A-coefficients. The position of the crossings depends also on the ac-Stark shift of the 
^5*0 level, which is dominated by the 423 nm ^5*0 — ^ Pi transition. But its Einstein coefficient is well known from 
photoassociation spectroscopy [24] and the shift of the ^5*0 level is therefore well known in this wavelength region. 
Unfortunately, the published Einstein coefficients of the 1.9 /zm line differ up to a factor of three (see Tab. I). 



III. EXPERIMENT 

To determine the light-induced shift of the ^5*0 — '^ Pi intercombination line in ^"Ca one has to compare its transition 
frequency in presence of the perturbing laser field to the frequency in the non-perturbed case. We measured the tran- 
sition frequency on a sample of laser-cooled atoms in a magneto-optical trap (MOT) by time-domain Ramsey-Borde 
atom interferometry. The experimental setup as well as the use of atom interferometry for the realization of an optical 
frequency standard were described previously [3, 4, 25-27]. 

In the experiment, lO'' calcium atoms were loaded from a thermal beam into a MOT and laser cooled within 20 ms 
to about 3 mK on the strong 423 nm transition from the ground state to the ^Pi level. The quadrupole field and 
laser beams of the MOT were then switched off, a homogeneous magnetic bias field was applied to separate the 
Zeeman components of the triplet state by 3.8 MHz each, and the intercombination line was interrogated by two 
counter-propagating laser pulse pairs during the expansion of the cold atomic cloud. The pulse separation within 
one pair was 216.4 fis. The light was generated by a diode laser in Liftman configuration, which was locked by the 
Pound-Drever-Hall method to a reference resonator. A line width of below 5 Hz and frequency drift below 0.5 Hz/s 
were achieved. This laser served as master laser, whereas for each direction of pulse pairs, an injection-locked slave 
diode was used to increase the power. The interrogation laser was stabilized to the central fringe of the Ramsey-Borde 
interference pattern (FWHM: 1.16 kHz) by adjusting the frequency offset between laser and reference cavity with an 
acousto-optical modulator. 

To measure the frequency shift due to the perturbing laser, we interleave two stabilization schemes with and without 
perturbing laser. The difference between the respective offset frequencies is the ac-Stark shift Av of the intercombina- 
tion line. This method suppresses systematic errors, e. g. frequency shifts due to misalignment of the interferometer. 
To investigate the dependence on the wavelength of the perturbation, radiation from different lasers was used. Near 
780 nm, we applied a diode laser with semiconductor amplifier (P = 300 mW at the atomic cloud). A multi-mode 
fiber-couple diode bar was used at 810 nm (P = 900 mW). At 860, 890, and 922 nm, a Ti:Sapphire laser delivered 
around P — 250 mW, and for 1064 nm, a Nd:YAG laser was available (P — 400 mW). The laser power was measured 
with wavelength-insensitive thermal power meters, which were calibrated at 800 nm. 

The radius w (1/e^ in irradiance) of these perturbing laser beams was typically in the order of 1 mm, while the 
Icr-radius R of the atomic cloud at the end of the trapping phase was about 0.4 mm. Under these conditions observed 
shifts are in the order of 10 Hz. The size of the beam is a compromise between high irradiance to achieve large 
frequency shifts and a moderate spatial inhomogeneity of the perturbation in the region of the atomic cloud. To 
take into account the residual inhomogeneity of the irradiance profile of the laser beam, we averaged over the beam 
and trap geometry. The laser beam was described as propagating in x-direction with constant radii Wy and Wz and 
Gaussian profile I{y, z) in radial direction. The measured density of the atomic cloud was well described as three 



dimensional Gaussian profile D{x, y, z) with R^ = Ry ^ Rz ■ The observed frequency shift Av is 

_ / dx dy dz s{u!, p, mp)I{y, z)D{x, y, z) 
J dx dy dz D{x, y, z) 

= s(uj,p,mp)-2P/{n{wl+4Rl)^/\wl+4Riy^'). (2) 

Here, we utilize the function 

s{uj,p,mp) = [ap{uj,p,mp) - as{LO,p,0)] /{2eohc), (3) 

which describes the frequency change of the intercombination line per unit irradiance of the perturbing field. 
Eq. 2 describes well the functional dependence on the geometrical quantities, as has been confirmed experimentally. 
However, it does not take into account the geometry of the spectroscopy and detection laser beams and the expansion 
of the atomic cloud. All influences were modeled by a Monte-Carlo simulation, leading to an additional correction 
depending on the actual conditions, but being always below ten percent. 

The experimental setup required to have the perturbing beam propagating perpendicularly to the quantization axis 
given by the magnetic bias field. In this geometry, linearly polarized light with its electric field vector parallel to the 
quantization axis is described as 7r-polarization, while with 90° rotated polarization one realizes cr-polarization, which 
can be decomposed into cr-polarization. 

For the ms = ^ mp = transition and 7r-light we find an easy to interpret situation, while for a two field 
components have to be considered simultaneously. For symmetry reasons of the 3 J-symbol in Eq. 1, a^- and cr^ -light 
generate the same shifts on the to = levels. Hence, the sum of these shifts is the same as for pure a^- or cr^-light. 
Fluctuations of the bias field prevented the direct observation of frequency shifts in the 10 Hz range of the mp = ±1 
Zeeman levels. Therefore, the spectroscopy laser was locked to the cross-over signal between the transitions to the ±1 
Zeeman components. For this purpose, the bias filed was reduced to allow for overlapping Doppler profiles of both 
Zeeman levels. Nevertheless, the measurements showed increased frequency noise. The measured frequency shifts 
for cr-polarized light determined in this way can be calculated as the average of the shifts for the transitions to the 
mp = 1 with (T+- and cr^-light. From measurements of the mp — 1 cross-over with 7r-light no new information would 
be gained here, because the ac-Stark shift is the same as for mp ~ with cr-light. 

IV. RESULTS AND DISCUSSION 

The experiments described in the previous section yielded the values of s{X,p,mp) depicted in Fig. 3. The uncer- 
tainties are due to noise of the measured frequency shift, uncertainty of the measured laser power (20 mW), uncertainty 
of the determination of the atomic cloud size R and beam radius w (total contribution of 10 %), and possibly non 
optimal overlap of laser beam and atomic cloud (10 %). 

To describe the observed ac-Stark shifts and determine the "magic wavelength" by a physical meaningful interpola- 
tion function (Eq. 3), the Einstein coefficients of the 1.9 /xm and 612 nm transitions were fitted to the measured data 
(Fig. 3). The former was chosen since it is not very well known (see Tab. I). It would have been sufficient to change 
this value to match the zero-crossing observed for the mp = 0, a case, but the absolute values of s, especially for 
mp — 0, TT would have been poorly described. Thus, A{612 nm) was also fitted to the data. As discussed at the end 
of Sec. II, the interplay of ^(612 nm) and ^4(1.9 /im) is crucial for the wavelength, at which the effective ac-Stark 
shift of the intercombination line vanishes. Other transitions are far detuned and lead to a weak dependence of the 
Stark shift on the wavelength in the investigated wavelength interval. Consequently, few information can be gained 
about other transitions. 

Our experimental data is very well reproduced by the fit. In order to estimate the uncertainties of the so obtained 
"magic wavelength" and Einstein coefficients a Monte-Carlo simulation approach was used. Artificial data sets were 
generated with a statistical distribution around the values calculated from the best fit and a width of the distribution 
according to the experimental uncertainty. An additional random parameter with average one and standard deviation 
cr = 0.1 was used as common scaling factor for all data to account for systematic uncertainties like the measurement 
of laser power. Then, these data sets were fitted by Eq. 3 and the obtained values for the A^i's and the "magic 
wavelength" were utilized for the error analysis. 

Fig. 4 shows the fitted Einstein coefficients of the fine structure multiplets A{1.9 /im) and A(612 nm) of 2000 such 
iterations. A clear correlation between both parameters is visible. The combined 68.3%-confidence limits for both 
parameters are also shown. The fitted Aki^s are listed in Tab. I together with values from literature. All uncertainties 
correspond to one standard deviation. Both A coefficients agree within the uncertainties with theoretical [15-18] 
and experimental [14, 19, 23] values from other authors. For the "magic wavelength" Aj ^ we found 800.8(22) nm. 



Subscripts denote mp and the polarization, the superscript indicates the angular momentum J of the ^Pj state. The 
uncertainty is one standard deviation of the distribution of values A calculated from the artificial Monte-Carlo data 
and parameter sets. 

Furthermore, the fitted parameters were cross checked by calculating the tensor polarizability atons of the ^Pi level 
and the difference Actscai of the scalar polarizabilities of this and the ^5*0 state, from the expressions given by Angel 
and Sanders [28]: 

^ 2TTeoc^2Jk + 1) 



/ 40J,(2J,-1) ^^ ^^j^_j^ \ 1 1 2 1 3^6oc^(2Jfc + 1^ ^ 
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The expression in curly brackets is a 6J-symbol, analytical expressions are tabulated in e.g. [11]. For both val- 
ues an uncertainty analysis as described above was performed. We found Aascai — 33(8) kHz/(kV/cm)^ and 
tttens = 2.9(6) kHz/(kV/cm)^. Within the uncertainties, they agree well with atens — 2.623(15) kHz/(kV/cm)^ 
as measured by Yanagimachi et al. [23] and the value of Aascai — 29.874(87) kHz/(kV/cm)^ one can calculate with 
atens from [23] and the Stark shift measurements of Li and van Wijngaarden [29]. Therefore, this check confirmed 
the results of our fit, but it is not suited to determine Aascai or atens more precisely. 

If these polarizabilities with their small uncertainties are included in the fit of the parameters yl(612 nm) and 
A(1.9 /im), these coefficients are almost uniquely determined by the static polarizabilities and the ac-Stark data 
is not well described. Therefore in principle the Einstein coefficients of all transitions should be fitted to the data, 
which however is not feasible, as it would lead to an underdetermined set of least-square equations. To estimate the 
influence of other transitions, we have included one of the Einstein coefficients of the strong transitions from the ^P 
level to the levels 4s4ii '^D (443 nm), Ap^ ^P (430 nm), icf ^P (300 nm), or Asbd ^D (363 nm) in the fitted parameters. 
E.g. including the transition to the 4s4d '^D (443 nm) level in the fit renders a good description of Aascai, atens, and 
of our measurements. ^(443 nm) has to be changed by only -4 % with respect to the value given in [14]. Similar 
results are obtained by fitting the other coefficients. The necessary variations are -1.5 %, -9 %, and -20 %, respectively. 
When we included more lines in the fit, no meaningful results could be obtained because of the least-square problem 
starts to become underdetermined. 

Values of Aq ^, A(612 nm), and ^(1.9 /im) determined this way do not depend on the choice of the additional line 
and do not differ within the uncertainties from the above discussed results. Only the uncertainty of the Einstein 
coefficients is reduced. This is due to a weaker correlation between ^(612 nm) and yl(1.9 /im) (compare Fig. 4), since 
A(1.9 /im) strongly influences atens- The average values of both coefficients obtained in such fits are listed in the 
second row of Tab. I. 

We have used the fitted A's to calculate the second "magic wavelength" of the "^Pi, A^j^ ^_ , and the one for the 
^Pq level, Aq, to which the transition is allowed in the fermionic isotope. Both were not observed directly in the 
experiment. The uncertainties were determined in the same way as for Aq^.. The determined values are given in 
Tab. II (uncertainties are one standard deviation), examples of the number distributions of the error analysis are 
shown in Fig. 5. The results of fits with or without including the polarizabilities in the dataset do not differ within 
the uncertainties. 

V. IMPLICATIONS FOR AN OPTICAL LATTICE CLOCK 

Because the light shift in the ^Pq state of the fermionic isotope '''^Ca is only very weakly dependent on the polarization 
of the lattice laser [8], it seems to be the best candidate for an optical lattice clock. However, the low natural abundance 
(0.135%) and its half-integer nuclear spin (/ — 5/2) that leads to magnetic-field sensitive transitions could become a 
problem. For the use of the "^Pi state in an optical lattice clock the transition is dependent on the polarization of the 
lattice laser with respect to the magnetic quantization field. For a dipole laser of power 1 W focussed to a waist of 
30 yum a trap depth of 35 /iK is achieved with cr-polarization at the "magic wavelength" at 800 nm. The shift for the 
"wrong" TT-polarization amounts to <5^ = 200 kHz. For a typical magnetic bias field of P = 3 /iT the magnetic field 
mostly determines the quantization axis, and the shift varies as (5^ sin {(p), where (f> is the angular deviation from the 
perpendicular direction. An experimentally achievable deviation of = 10 mrad would amount to a shift of 20 Hz. 
However, as this shift depends on the power of the trap laser, it can be identified and corrected for, e. g. by a small 
change in the wavelength of the trapping laser. We have calculated the shift for the combined electric and magnetic 
fields by diagonalizing the total Hamiltonian. For the example given above we find that at a frequency shift of 7 GHz 
of the trapping laser the derivative of the light shift on the trapping laser power vanishes. Varying the power also 



changes the effective quantization axis for the combined E and B-fields. Thus the derivative depends on the power and 
the Hght shift at this point is not precisely canceled. But even for a relatively small field of 3 /zT the quantization axis 
is mainly defined by the magnetic field and the resulting frequency offset from the unperturbed line is only 40 mHz. 
Similar shifts have been found for residual ellipticity of the trapping laser beam of 4 x 10^'"'. Therefore also the ■^Pi 
state could be employed in an optical lattice clock where the residual AC-Stark shift contributes less than 10^^^ to 
the uncertainty. 

VI. CONCLUSION 

We have measured the ac-Stark shift of the ^"Ca intercombination line and analyzed our data by modeling the 
wavelength dependence. Our data is well described after fitting the Einstein coefRcients A of the transitions 4s3(i ^D — 
AsAp ^P and 4s5s '^5' — AsAp ^P to our data (Fig. 3). The A coefficients agree with experimental results [14, 19, 23] 
and confirm recent theoretical calculations [15-18] of the oscillator strengths (Tab. I). The analysis allowed for the 
determination of wavelengths, at which the energy shifts of both levels of the intercombination line are equal, the so 
called "magic wavelengths" (Tab. II). We found "magic wavelengths" at 800.8(22) nm and 983(12) nm for the ^Pi 
level and 735.5(20) nm for the transition to the ^Po level, which are of interest for optical lattice clocks. We have 
proposed a method that is capable to reduce the influence of the trapping field to below 10~^^ even with moderate 
requirements on the quality of its polarization. 
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FIG. 1: Simplified level scheme of the lowest levels in *'^Ca and couplings by the perturbation laser. The energies are not to 
scale. 
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FIG. 2: ac-Stark shift V for a laser irradiance of 1 W/m^. Some combinations of Zeeman levels and polarizations show the 
same shift. The poles at 612 nm in two curves of the ^Pi manifold and the ^Po curve are due to coupling to the 4s5s ''5' level. 
Circles denote points of vanishing frequency shift of the intercombination line, so called "magic wavelengths". 
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FIG. 3: Measured frequency shift function s measured as described in Sec. Ill (symbols). Lines are the best fit according to 
the model introduced in Sec. II. For details about the fitting procedure see Sec. IV. 
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FIG. 4: Einstein coefficients A{^D — ^P) and A{^S — ^P) fitted to 2000 artificial data sets generated by a Monte-Carlo approach. 
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FIG. 5: Number distribution for the "magic wavelengths" \q^^ and A^-^ ^_ as determined by the Monte-Carlo simulation. 



TABLE I: Determined Einstein A coefficients (first two fines) and values from literature. The last column indicates theoretical 
or experimental origin of data. 



reference 



A{4sM'-^D - 4s4p^P) 



A{4s5s'-^S - 4s4p'^P) 



origin 



this work " 
this work ' 

[12] 
[14] 
[15] 
[16] 
[17] 
[18] 
[19] 
[23] 



8.6(25)-10'' s 
7.8(4)-10'^ s" 

s" 



g.i-io-' s"^ 
s-io""' s-i 

8.4(25) •lO-'"' s^i 
7.06 -lO' s^i 



8.7(8)-10' s" 
8.5(4)-10^ s 



6.6 

8.6 
8.3 
8.1 

7.8 



10' s" 
10^ s" 
10'^ s" 



10^ 



8.6(9)-10^ 
8.6(13)-10^ 



exp. 
exp. 
theo./exp. 
exp. 
theo. 
theo. 
theo. 
theo. 
exp. 
exp. 



"Fitted without static polarizabilities AcVscali citens- 

''Fitted including static polarizabilities and one additional line (see text). 



TABLE II: "Magic wavelengths" 
belongs to the perturbing laser. 



determined in this work, mp denotes the upper state's Zeeman level, the polarization given 



transition 



rnp 



4s^ "So - 4s4p ^Pi 
4s2 ^So - 4s4p ^Pi 
4s2 ^So - 4s4p ^Po 



polarization 



value 








a- (a+) 
arbitrary 



800.8(22) nm 

983(12) nm 

735.5(20) nm 



